Stochasticity effects in quantum radiation reaction 
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When an ultrarelativistic electron beam collides with a sufficiently intense laser pulse, radiation- 
reaction effects can strongly alter the beam dynamics. In the realm of classical electrodynamics, 
radiation reaction has a beneficial effect on the electron beam as it tends to reduce its energy spread. 
Here, we show that when quantum effects become important, radiation reaction induces the opposite 
effect, i.e., the electron beam spreads out after interacting with the laser pulse. We identify the 
physical origin of this opposite tendency in the intrinsic stochasticity of photon emission, which 
becomes substantial in the full quantum regime. Our numerical simulations indicated that the 
predicted effects of the stochasticity can be measured already with presently available lasers and 
electron accelerators. 
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A deep understanding of the dynamics of electric 
charges driven by electromagnetic fields is one of the most 
fundamental problems in physics, as it has implications 
in different fields, including accelerator, radiation and 
high-energy physics. Apart from its impact on practical 
issues, as the construction of new experimental devices 
(e.g., quantum x-free electron lasers [1]), the investiga- 
tion of the dynamics of electric charges (electrons, for 
definiteness) is also of pure theoretical interest, as it in- 
volves in general a coupled dynamics of the electrons and 
of their own electromagnetic field. 

In the realm of classical electrodynamics, radiation- 
reaction (RR) effects stem from the back reaction on 
the electron dynamics of the electromagnetic field gen- 
erated by the electron itself while being accelerated by 
a background electromagnetic field [3, 'Sj . The Landau- 
Lifshitz (LL) equation has been recently identified as the 
classical equation of motion of an electron, with mass 
m and charge e < 0, which includes RR effects self- 
consistently [2- 7] , although alternative models have been 
suggested [8|, IH . The analytical solution of the LL equa- 
tion in a plane-wave field [10] shows that if an electron 
impinges with initial four-momentum onto a plane- 
wave field (electric-field amplitude Eq^ central angular 
frequency uoq and propagating along the direction n), 
RR effects substantially affect the electron dynamics, 
if the parameter Rc = Q^Xo^o is of the order of unity 
(see also [H]). Here, a = is the fine-structure con- 
stant, xo = {{npo)/Tn)Eo/Ecr is the so-called quantum 
nonlinearity parameter, with = (l,n) and Ecr = 
mV|e| = L3 x iQiev/cm^ ^ \e\Eo/muJo is the 

classical nonlinearity or relativistic parameter (units with 
h = c = 1 are used throughout). It is worth noting that, 
although Xo is much smaller than unity in the realm of 
classical electrodynamics [2L the parameter Rc can be 
of the order of unity [J, M lUI]- 

The parameter Rc 
represents the average energy radiated by the electron 
in one laser period in units of the initial electron en- 
ergy, and for an ultrarelativistic electron initially coun- 
terpropagating with respect to the laser field with en- 



ergy £, it is Rc = 3.2£:[GeV]/o[10^^ W/cm2]/cjo[eV], with 
Iq = Eq/ Att being the laser pulse peak intensity. The nu- 
merical value of the parameter Rc shows the generally 
demanding requirements to observe large RR effects and 
it explains why the LL equation still lacks an experi- 
mental confirmation (see 11- 13] for recent experimental 
proposals). The expression of the parameter Rc is also in 
agreement with the well-known classical result that more 
energetic particles radiate more at given other conditions 
[15]. In turn, this explains physically the beneficial effect 
of RR when it is included, e.g., in the investigation of 
the production of electron [16] and ion [rzLfioj bunches 
in laser-plasma interaction. In fact, it is found that RR 
acts as a cooling mechanism and its effects render the 
energy spectra of the produced particle bunches more 
monochromatic than if RR is not included. 

In this Letter we show that when quantum effects be- 
come important RR induces exactly the opposite behav- 
ior and makes the energy distribution of an electron beam 
initially counterpropagating with respect to a strong laser 
field broader as it was before the interaction. We explain 
this striking difference between classical and quantum 
RR relating it to the stochastic nature of the emission 
of radiation, which becomes substantial in the quantum 
regime, and indicating that quantum effects amount to 
add a stochastic term in the LL classical equation. By 
means of numerical simulations we show that the broad- 
ening of the electron distribution in the quantum regime, 
is measurable in principle with presently available tech- 
nology even in an all-optical setup. Our results are rele- 
vant for future laser-based electron accelerators, indicat- 
ing that one cannot rely on the beneficial effects of RR 
on the energy spread of the electron beam at sufficiently 
high electron energies that quantum effects become im- 
portant. We note that the stochastic nature of photon 
emission has instead been shown to lower the laser in- 
tensity threshold at which electromagnetic cascades are 
generated [2l|. 

Taking into account exactly RR in the full strong-field 
QED regime is a formidable task, as it amounts to deter- 
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mine completely the S'-matrix, describing the interaction 
of the electron-positron field with the radiation field in 
the presence of the strong background electromagnetic 
field [2|, [22|. Thus, we limit here to the so-called "nonlin- 
ear moderately-quantum" regime fiij, where: 1) ^ 1, 
such that nonlinear effects in the laser field amplitude 
are large; 2) xo ^ 1, such that nonlinear QED effects 
are already important, but electron-positron pair pro- 
duction is still negligible. In this regime, RR effects on 
the electron dynamics in a strong plane- wave field mainly 
stem from the sequential emission of many photons by 
the electron, and they can be investigated by means of 
a kinetic approach [23l425j (see [22|, for an alternative, 
microscopic approach). In this approach, the electrons 
and the photons are described by distribution functions 
in phase space, which obey to "kinetic" equations. Since 
electron-positron pair production is neglected: 1) the dis- 
tribution function of positrons can be assumed to vanish 
identically; 2) the kinetic equation for the electron dis- 
tribution function is not coupled to that of the photons 
(isl-liil. Another realistic approximation, which allows 
us to avoid technical complications in favor of a clearer 
physical understanding, is to consider an electron bunch 
initially counterpropagating with respect to the laser field 
and with a typical energy e* m^o- This is the case, 
for example, in the realistic situation of a electron bunch 
with typical energy e* = 1 GeV colliding head-on with an 
optical {ujQ = 1.55 eV) laser field of intensity 10^^ W/cm^ 
[26] for which m^o = 25 MeV. The condition e* > m^o 
ensures that the transverse momentum of the electrons 
(with respect to the initial propagation direction) re- 
mains much smaller than the longitudinal one in passing 
through the plane wave [13] , and this reduces the present 
problem to a one-dimensional one. 

By assuming that the plane wave propagates along 
the positive y direction and that it is linearly-polarized 
along the z direction, we can write its electric field as 
E{(p) = Eof{(p)z^ where cp = uJo{t — y) is the laser 
phase and f{(p) is the pulse-shape function such that 
1/(^)1 max = 1- If = {s^p) is the four- momentum 
of an electron, it is convenient to introduce the quantity 
P- = € — Py^ which is a constant of motion in the plane- 
wave field under consideration [27^. However, if the elec- 
tron emits a photon with four-momentum = (oJ^k)^ 
then its four-momentum changes to p^^ = {e^p') and 
y_ = p- — with p^_ = e' — p'y and k- = u — ky. The 
single-photon emission probability per unit phase (p and 
per unit u = k-/{p- — k-) in the ultrarelativistic regime 
^0 > 1 reads [28| 



1 



f 2u 



1 



V3x(^,P- 



/oo 
dxY^i {x) 
2u 3 
■. ^,/ 7 



3x(¥',P_ ) 



where Kjy(-) is the modified Bessel function of z/th order 
and where x{^iP-) = (P-/^)|^(v^)|/^cr, with E{(p) = 
Eof{(p). Since the probability in Eq. ([1]) depends only 
on the phase-space variables cp and p_, it is possible to 
describe the electron beam via an electron distribution 
ne((f^P-)^ which satisfies the kinetic equation (see Ref. 
[24]) 
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The kinetic equation (j2j) is an integro-differential equa- 
tion, i.e., it is non-local in the momentum p-. This oc- 
currence is intimately connected to the quantum nature 
of the emission of radiation. In fact, the emission of ra- 
diation is described quantum mechanically as the emis- 
sion of photons, which carry energy and momentum, such 
that, if an electron emits a photon with momentum 
its initial state with a given momentum po,- will be cou- 
pled to that with momentum po,- — k-^ with k- ranging 
from to po,-- 

In order to investigate the classical limit of Eq. (j2j) 
for xi^^P-) <^ 1? it is convenient to perform the change 
of variable v = {p'_ — p_)/p_x((^,p_) {v = k-/{p- — 
k-)x{^-,P-)) in the first (second) integral in Eq. ([2]). 
By expanding the resulting equation in x{^->P-) ^nd by 
keeping terms up to the order one obtains the 

Fokker-Planck-like equation jUj (see also (25I. [3Q|) 
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with a "drift" coefficient A{(p^p_) and a "diffusion" co- 
efficient B{(p^p-) given by 
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respectively. It is worth observing that this equation is 
no longer an integro-differential equation but rather a 
partial differential equation. In other words, when quan- 
tum photon-recoil effects become smaller and smaller, the 
distribution function of electrons with a momentum p- 
depends only on its values close to p- and its dynamics 
is local. On this respect, we also note that higher-order 
corrections in xi^^P-) would result in the appearance 
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of terms proportional to higher and higher derivatives of 
ne{(p,P-) with respect to p_. 

If we first consider only the terms proportional to 
X^i^^P-) in Eq. ([5]), the latter equation has the form 
of a Liouville equation: 



dip 



d 
dp- 



dp- 



where 



dp- 

dif 



Iclj^.P- 



(8) 



(9) 



with Ici{^^P-) = (2/3)am^y^((p, p_) being the classi- 
cal intensity of radiation [27|. Equation ([9]) is exactly 
the classical single-particle equation for the momentum 
P- resulting from the LL equation [10] (see also [31]). 
In other words, the terms in Eq. (j5j) proportional to 
X^{^^P-) describe the classical dynamics of the electron 
distribution including RR. The fact that Eq. (|8]) has 
the form of a Liouville equation implies, as it must be, 
that the classical dynamics of the electron distribution 
is deterministic [29|]. Also, since the single-particle equa- 
tion (j9j) admits the analytical solution [10.], p-{(p;po^-) = 
Po-/h{(p,po-) with 
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for an electron with initial momentum p^(0) 



= = 



(eo,Po) (Po,- = Co — p^^y\ one can write explicitly the 
exact analytical solution of Eq. (|8]) by means of the 
method of characteristics. If the distribution ne(0,p_) at 
the initial phase = is given, for example, by the Gaus- 
sian distribution ne(0,p_) = A^exp[— (p_ — P^)^/2cr^_], 
where is a normalization factor, pi is the average value 
oip- and crp_ is the standard deviation [32], then the so- 
lution of Eq. (jH]) reads 
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with g{(p^p_) = h{(p^ —p-). Since po,- in Eq. (p!Q|) is 
positive for finite values of po^y and po,- ^ only at 
Py +00, the function g{(p^p_) must be non- negative 
for all values of and the equation g{(p^P-^max) = 
fixes the maximum value p-,max = P-,max(v^) allowed 
for the variable p- at each (p. Before passing to in- 
vestigate the quantum corrections in Eq. (|5]), we ob- 
serve that the classical solution in Eq. (p!Q|) is such 
that < dp-{(p; po^-) / dpo^- < 1 for > and this 
ensures that, due to RR effects, the difference Ap-{(p) 
between the momenta of two electrons decreases for in- 
creasing values of (p. This implies that RR effects tend 
to decrease the energy width of the electron distribu- 
tion in agreement with previous results 19|, [20|. Also, if 



o'p_ <C pl_ in Eq. (pTj) . it can be seen that the distribu- 
tion ne{(p,P-) is approximately a Gaussian with effective 
width ap_{(p) ^ ap_ /h'^{(p^p*_) decreasing at increasing 

The quantum corrections in Eq. ^ to the classical 
kinetic equation (|8]) stem from two different contribu- 
tions. The first one affects the drift coefficient A{(p^p-) 
(see Eq. ([6])) and it corresponds to the leading quan- 
tum correction to the total intensity of radiation found in 
0, [m . This correction, does not change the structure of 
the classical kinetic equation (|8]) but only the "effective" 
momentum change per unit phase. Since this leading 
quantum correction is negative, we expect that it tends 
to decrease the reduction of the width with respect to the 
classical prediction. However, by replacing the classical 
intensity of radiation Ici{(p^P-) with the corresponding 
quantum one Iq{(p^p-) (see, e.g., Eq. (83) on pg. 522 
in [28]), the resulting Liouville equation would still pre- 
dict a reduction of the width of the electron distribution 
function. Although this statement can be proven math- 
ematically, it can be intuitively understood as a physi- 
cal consequence of the fact that more energetic electrons 
on average emit more radiation. On the other hand, 
however, the second leading quantum correction corre- 
sponds to the diffusion coefficient B{(p,p-) in Eq. ([7]) and 
it alters the structure of the classical kinetic equation. 
The appearance of a diffusion- like term in the kinetic 
equation of the electron distribution is intimately con- 
nected to the stochastic nature of the quantum emission 
of photons. According to the theory of stochastic differ- 
ential equations, in fact, the Fokker- Planck- like equation 
is related to the single-particle stochastic equation 
dp- = —A{(p^p-)d(p + B{(p^p-)dW ^ where dW repre- 
sents an infinitesimal stochastic function [29]. The diffu- 
sion term in Eq. (j5j) is responsible of the broadening of 
the distribution function. In the case, for example, of a 
Gaussian distribution function assumed to be well peaked 
at (p around the classical value p-{(p]p*_) (see Eq. (fTO]) ). 
it can easily be shown, that if cFp_ is its initial width, 
then 
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This result clearly shows the opposite influence of the 
classical drift term and of the quantum diffusion term on 
the width of the electron distribution. It is worth noting 
that the correction to the width arising from the quan- 
tum corrections in the drift coefficient A{(p^p) is found 
to be smaller than the correction proportional to the dif- 
fusion term by a factor cr'^/p*!^ <C 1, and it has been 
neglected in the approximated expression ([12]) . We warn 
the reader about the fact that a formal solution of the 
Fokker-Planck equation, for example, in the case of an 
initial ^-like momentum distribution and vanishing drift 
term, would predict, due to the spreading in the mo- 
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mentum distribution, the spurious presence of particles 
with momentum larger than the initial one. This indi- 
cates that, for a completely consistent treatment, the full 
equation ([2]) has to be employed, which will be carried 
out below numerically. 

The above effect on the broadening of the electron mo- 
mentum distribution can also be interpreted in terms of 
the entropy S{(f) = - dp- ne{(f,P-) \ii[ne{(f,P-)/no] 
associated to the the electron distribution, where the 
Boltzmann constant has been set equal to unity and 
where the physically ineffective constant no can be cho- 
sen, for example, such that S{0) = 0. By employing this 
definition and Eqs. (|5]), (|6]) and (O, it results 
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This result further corroborates the idea that, while the 
classical "deterministic" evolution of the electrons im- 
plies a reduction of the entropy of the electrons, quantum 
corrections tend to increase it. 

In order to show that the effects discussed above can be 
in principle measured with presently available laser and 
electron accelerator technology, we consider below two 
numerical examples. In both cases we assume a laser 
pulse with f{(f) = sm'^{(p/2NL) sm{(p) for cp G [0, = 
[0,2A^L7r] and zero elsewhere, where Nl is the number 
of laser cycles and with ujo = 1-55 eV, and an initial 
Gaussian electron distribution with a total number of 
1000 electrons. 

In the first numerical example, we choose the laser and 
electron parameters such that quantum effects are neg- 
ligible, whereas RR effects are relatively large. We set 
Jo = 4.3 X 10^0 W/cm^, p^_ = 84 MeV (note that £* ^ 
pl/2 = 42 MeV) such that x* = {pl/m){Eo/Ecr) - 
5 X 10"^, ap_ = 8.4 MeV, and Nl = 1600, correspond- 
ing to a pulse duration of about 4 ps. The results for 
the initial and final distribution are shown in Fig. 1. As 
expected, the final distribution ne{(pf^p-)^ calculated by 
solving numerically Eq. (|2]) (solid, red line) and the clas- 
sical analytical solution n^^{(pf,p-) (see Eq. (pT]) ) are 
very similar and both show a reduction of the width from 
the initial value 8.4 MeV to the final one 4.7 MeV. In 
the second numerical example, instead, we want to probe 
the quantum regime and we set /q = 2 x 10^^ W/cm^ 
[26], pi = 2 GeV (e* ^ 1 GeV) and ap_ = 0.2 GeV 
corresponding to x* = O-S^ ^ind Nl = 10 correspond- 
ing to about 30 fs. Electron beams with such energies 
are nowadays available not only in conventional acceler- 
ators but also by employing plasma-based electron accel- 
erators Q (see also ^]), allowing in principle for an 
all-optical setup. The results of our numerical simula- 
tions are shown in Fig. 2. The figure shows that in the 
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FIG. 1. Comparison of the initial electron distribution (dot- 
ted, blue line) and the final electron distribution according to 
Eq. ([2]) (solid, red line) and to Eq. (pT]) (dashed, green line). 
The laser and the electron distribution parameters are given 
in the text. 
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FIG. 2. Time evolution of the electron distribution for a 10- 
cycle sin^-like laser pulse (part a)) according to Eq. ([2|) (part 
b)), to Eq. {[n]) (part c)) and to Eq. JH)) and Eq. ^ with the 
replacement Ici{^,P-) Iq{(p,p-) (part d)). The laser and 
the electron distribution parameters are given in the text. 



quantum regime the full quantum calculations based on 
Eq. ([2]) predict a broadening of the electron distribution 
(Fig. 2b), according to our analysis above. Whereas, 
the classical calculations based on the exact solution in 
Eq. (fTTj) (see Fig. 2c) predict a strong narrowing of the 
distribution. It is interesting to note that, according to 
the discussion above Eq. (p!2]) . if we consider the classical 
equation ([8j) and we substitute the classical intensity of 
radiation /^/(y, p ) with the quantum intensity Iq{(p^p-) 
(see, e.g., [24, 28]), the corresponding results (see Fig. 
2d) still predict a narrowing of the distribution function. 
This clearly supports the idea that the broadening of the 
electron distribution is an effect of the importance of the 
stochasticity of the emission of radiation, which becomes 
substantial in the quantum regime. 

In conclusion, we have demonstrated that the impor- 
tance of the stochastic nature of the emission of radia- 
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tion in the quantum regime, has a profound impact on the 
evolution of an electron beam passing through an intense 
laser field. The stochasticity, in fact, induces a broad- 
ening on the electron momentum distribution, whereas 
classical theory of RR even predicts a narrowing of the 
distribution itself. A numerical example has shown the 
feasibility of measuring such effects by employing already 
demonstrated laser intensities and electron-beam ener- 
gies. 

The authors would like to gratefully acknowledge use- 
ful discussions with John Kirk and Tom Blackburn. 
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